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Abstract 

Inclusive transverse momentum spectra of primary charged particles in Pb-Pb collisions at ^/s^ 
= 2.76 TeV have been measured by the ALICE Collaboration at the LHC. The data are presented 
for central and peripheral collisions, corresponding to 0-5% and 70-80% of the hadronic Pb-Pb 
cross section. The measured charged particle spectra in [tjI < 0.8 and 0.3 < pr < 20 GeV/c are 
compared to the expectation in pp collisions at the same y/s^, scaled by the number of underlying 
nucleon-nucleon collisions. The comparison is expressed in terms of the nuclear modification factor 
Raa- The result indicates only weak medium effects (Raa ~ 0.7) in peripheral collisions. In central 
collisions, Raa reaches a minimum of about 0.14 at pr = 6-7 GeV/c and increases significantly at 
larger pr- The measured suppression of high-/??' particles is stronger than that observed at lower 
collision energies, indicating that a very dense medium is formed in central Pb-Pb collisions at the 
LHC. 



*See Appendix lAlfor the list of collaboration members 
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High energy heavy-ion coUisions enable the study of strongly interacting matter under extreme condi- 
tions. At sufficiently high collision energies Quantum-Chromodynamics (QCD) predicts that hot and 
dense deconfined matter, commonly referred to as the Quark-Gluon Plasma (QGP), is formed. With the 
advent of a new generation of experiments at the CERN Large Hadron Collider (LHC) Q a new energy 
domain is accessible to study the properties of this state. 

Previous experiments at the Relativistic Heavy Ion Collider (RHIC) reported that hadron production 
at high transverse momentum (px) in central (head-on) Au-Au collisions at a centre-of-mass energy 
per nucleon pair -^/S^ of 200 GeV is suppressed by a factor 4-5 compared to expectations from an 
independent superposition of nucleon-nucleon (NN) collisions E [3l IH 13. The dominant production 
mechanism for high-pr hadrons is the fragmentation of high-;??- partons that originate in hard scatterings 
in the early stage of the nuclear collision. The observed suppression at RHIC is generally attributed to 
energy loss of the partons as they propagate thi^ough the hot and dense QCD medium 1161171 [8ll9l[T0l. 

To quantify nucleai^ medium effects at high px, the so called nuclear modification factor Raa is used. 
Raa is defined as the ratio of the charged particle yield in Pb-Pb to that in pp, scaled by the number of 
binary nucleon-nucleon collisions (A'^coii) 

„ , , _ {l/N^)d^N^/d7ldpT 

(^coil) {l/N!^^)d^N^^'/d7ldpT ' 

where rj = — ln(tan0/2) is the pseudo-rapidity and 6 is the polai^ angle between the charged particle 
direction and the beam axis. The number of binary nucleon-nucleon collisions (A'coii) is given by the 
product of the nuclear overlap function (Taa) ifTTl and the inelastic NN cross section CJ;^^. If no nuclear 
modification is present, Raa is unity at high pj. 

At the larger LHC energy the density of the medium is expected to be higher than at RHIC, leading to a 
larger energy loss of high pj paitons. On the other hand, the less steeply falling spectrum at the higher 
energy will lead to a smaller suppression in the pr spectrum of charged particles, for a given magnitude 
of partonic energy loss ||9l[T0|. Both the value of Raa in central collisions as well as its pr dependence 
may also in part be influenced by gluon shadowing and saturation effects, which in general decrease with 
increasing x and Q^. 

This Letter reports the measurement of the inclusive primary charged particle transverse momentum 
distributions at mid-rapidity in central and peripheral Pb-Pb collisions at .^s^ = 2.76 TeV by the ALICE 
experiment |[T2ll . Primary particles are defined as prompt particles produced in the collision, including 
decay products, except those from weak decays of strange particles. The data were collected in the first 
heavy-ion collision period at the LHC. A detailed description of the experiment can be found in |[T2ll . 

For the present analysis, charged particle tracking utilizes the Inner Tracking System (ITS) and the Time 
Projection Chamber (TPC) fil3il . both of which cover the central region in the pseudo-rapidity range 
|t] I < 0.9. The ITS and TPC detectors are located in the ALICE central ban^el and operate in the 0.5 T 
magnetic field of a large solenoidal magnet. The TPC is a cylindrical drift detector with two readout 
planes on the endcaps. The active volume covers 85 < r < 247 cm and —250 < z < 250 cm in the radial 
and longitudinal directions, respectively. A high voltage membrane at z = divides the active volume 
into two halves and provides the electric drift field of 400 V/cm, resulting in a maximum drift time of 
94 ixs. 

The ITS is used for charged particle tracking and trigger purposes. It is composed of six cylindrical layers 
of high resolution silicon tracking detectors with radial distances to the beam line from 3.9 to 43 cm. The 
two innermost layers are the Silicon Pixel Detectors (SPD) with a total of 9.8 million pixels, read out by 
1200 chips. Each chip provides a fast signal if at least one of its pixels is hit. The signals from the 1200 
chips are combined in a programmable logic unit which supplies a trigger signal. The SPD contributes 
to the minimum-bias trigger, if hits are detected on at least two chips on the outer layer. The SPD is 
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Table 1: The average numbers of participating nucleons (A^pait), binary nucleon-nucleon collisions (A^coii)- and the 
average nuclear overlap function {Taa) for the two centrality bins, expressed in percentages of the hadronic cross 
section. 



Centrality 


(^part) 


(^coll) 


{TAA){mh-') 


0-5% 


383 ±2 


1690 ±131 


26.4 ±0.5 


70-80% 


15.4±0.4 


15.7 ±0.7 


0.25 ±0.01 



followed by two layers of Silicon Drift Detectors (SDD) with 133k readout channels. The two outermost 
layers are Silicon Strip Detectors (SSD) consisting of double-sided silicon micro-strip sensors, for a total 
of 2.6 million readout channels. 

The two forward scintillator hodoscopes (VZERO-A and VZERO-C) cover the pseudo-rapidity ranges 
2.8 < T] < 5. 1 and — 3.7 < T] < — 1 .7. The sum of the amplitudes of the signals in the VZERO scintillators 
is used as a measure for the event centrality. The VZERO detectors also provide a fast trigger signal if at 
least one particle hit was detected. 

During the heavy-ion data-taking period, up to 1 14 bunches, each containing about 7 x 10^ ions of ^''^Pb, 
were collided at .JT^ = 2.76 TeV in the ALICE interaction region. The rate of hadronic events was 
about 100 Hz, coiTcsponding to an estimated luminosity of 1.3x10^^ cm^^s^'. The detector readout 
was triggered by the LHC bunch-crossing signal and a minimum-bias interaction trigger based on trig- 
ger signals from VZERO-A, VZERO-C, and SPD. The present analysis combines runs taken with two 
different minimum-bias conditions. In the first set of runs, two out of the three trigger signals were re- 
quired, while in the second set a coincidence between VZERO-A and VZERO-C was used. Both trigger 
conditions have similar efficiency for hadronic interactions, but the latter suppresses a large fraction of 
electromagnetic reactions. 

The following analysis is based on 2.3 x 10^ minimum-bias Pb-Pb events, which passed the offline event 
selection. This selection is based on VZERO timing information and the correlation between TPC tracks 
and hits in the SPD to reject background events coming from parasitic beam interactions. Additionally, 
a minimal energy deposit in the Zero Degree Calorimeters (ZDC) is required to further suppress electro- 
magnetic interactions. Only events with reconstructed vertex at |zvtx| < 10 cm were used. The definition 
of the event centrality is based on the sum of the amplitudes measured in the VZERO detectors as de- 
scribed in |[T4l . Alternative centrality measures utilize the cluster multiplicity in the outer layer of the 
SPD or the multiplicity of reconstructed tracks. The correlation between the VZERO amplitude and the 
uncorrected TPC track multiplicity in [tjI < 0.8 is illustrated in Fig[T] The VZERO amplitude distribu- 
tion is fitted using a Glauber model ifTTIl to determine percentage intervals of the hadronic cross section, 
as described in llT4l . We used a Glauber model Monte Carlo simulation assuming aj^^ = 64 mb, a 
Woods-Saxon nuclear density with radius 6.62 ±0.06 fm and surface diffuseness 0.546 ±0.010 fm |[T5l . 
A minimum inter-nucleon distance of 0.4 ± 0.4 fm is assumed. The Glauber Monte Carlo allows one 
to relate the event classes to the mean numbers of participating nucleons (Npan) and binary collisions 
(A'^coii) (see Table [T]) by geometrically ordering events according to the impact parameter distribution. 
The eiTors include the experimental uncertainties in the parameters used in the Glauber simulation and 
an uncertainty of ±5 mb in CJj^^. The TPC multiplicity distributions for the central and peripheral event 
samples selected for this analysis, corresponding to the 0-5% and 70-80% most central fraction of the 
hadronic Pb-Pb cross section, are shown in the lower panel of Fig. [T] Charged particle tracks are recon- 
structed in the ITS and TPC detectors. Track candidates in the TPC are selected in the pseudo-rapidity 
range |t]| < 0.8. Track quality cuts in the TPC are based on the number of reconstructed space points 
(at least 70 out of a maximum of 159) and the per space point of the momentum fit (lower than 4). 
The TPC track candidates are projected to the ITS and used for further analysis, if at least two matching 
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Figure 1: Upper panel: Correlation between VZERO amplitude and the uncorrected track multiplicity in the TPC. 
Indicated are the cuts for the centrality ranges used in this analysis. Lower panel; Minimum-bias distribution of 
the TPC track multiplicity. The central (0-5%) and peripheral (70-80%) event subsamples used for this analysis 
are shown as grey histograms. 

hits in the ITS are found, including at least one in the SPD. The average number of associated hits in 
the ITS is 4.7 for the selected tracks. The event vertex is reconstructed by extrapolating the particle 
tracks to the interaction region. The event vertex reconstruction is fully efficient in both the peripheral 
and the central event sample. Tracks are rejected from the final sample if their distance of closest ap- 
proach to the reconstructed vertex in longitudinal and radial direction, and dxy, satisfies > 2 cm or 
<i.„ > 0.018cm + 0.035cm-;7j:^ °\ with pr in GeV/c. 

The efficiency and purity of primary charged particles using these cuts are estimated using a Monte Carlo 
simulation including HIJING |[T6l events and a GEANT3 lITTl model of the detector response [fT8). We 
used a HIJING tune which reproduces approximately the measured charged particle density in central 
collisions |[T4ll . In central events, the overall primary charged particle efficiency in |t]| < 0.8 is 60% at 
Pt = 0.3 GeV/c and increases to 65% at pj = 0.6 GeV/c and above. In peripheral events, the efficiency 
is larger by about 2-3%. The contamination from secondaries is 6% at pj = 0.3 GeV/c and decreases to 
about 2% at Pt > I GeV/c, with no significant centrality dependence. This contribution was estimated 
using the dxy distributions of data and HIJING and is consistent with a first estimate of the strangeness 
to charged particle ratio from the reconstruction of K^, A and A invariant mass peaks. 

The momentum of charged particles is reconstructed from the track curvature measured in the ITS and 
TPC. The momentum resolution can be parametrized as {o{pt)/ PtY = + {b- ptY- It is estimated 
from the track residuals to the momentum fit and verified by cosmic muon events and the width of the 
invariant mass peaks of A, A and K^. While a = 0.01 for both centrality bins, there is a weak centrality 
dependence of b, i.e. b = 0.0045 (GeV/c)^' in peripheral events and b = 0.0056 (GeV/c)^' in central 
events. This is related to a slight decrease for more central events of the average number of space points 
in the TPC. The modification of the spectra arising from the finite momentum resolution is estimated by 
Monte Cai4o. It results in an overestimate of the yield by up to 8% at pj = 20 GeV/c in central events. 
This was accounted for by introducing a pj- dependent correction factor to the pT spectra. From the 
mass difference between A and A and the ratio of positive over negative charged tracks, assuming charge 
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Table 2: Contributions to the systematic uncertainties on the inclusive spectra. For the pj dependent errors the 
ranges are given. 



Centrality class 


0-5% 


70-80% 


Centrality selection 


1% 


7% 


Track and event selection cuts 


1^% 


1^% 


Particle composition 




1^% 


Material budget 


1-2% 


1-2% 


Secondary particle rejection 


<1% 


<1% 


Tracking efficiency 


2-6% 


2-6% 


Total systematic uncertainties 


5-7% 


8-10% 



symmetry at high pj, the upper limit of the systematic uncertainty of the momentum scale is estimated 
to be \^{pt) / Pt\ < 0.002. This has negligible effect on the measured spectra. 

Table |2] shows the systematic uncertainties obtained by a comparison of different centrality measures 
(using the SPD instead of VZERO), and by varying the track and event quality cuts and the Monte Carlo 
assumptions. In particular, we studied a variation of the most abundant charged particle species (p, n, K) 
by ±30%, the material budget by ±7%, and the secondary yield from strangeness decays in the Monte 
Carlo by ±30%. We have used the differences between the standai^d analysis and one based only on 
the use of TPC tracks to estimate the uncertainty on the track efficiency corrections, included in the 
systematic errors. The total systematic uncertainties on the corrected pj spectra depend on pT and are 

8- 10% and 5-7% for the peripheral and central event samples, respectively. 

The determination of Raa requires a pp reference at y/s = 2.76 TeV, where no pp measurement exists. 
Different approaches are at hand which allow a prediction of the pj spectrum at a given y/s by scaling 
existing data at different energies. Such approaches assume general scaling properties of perturbative 
QCD (pQCD) or rely on next-to-leading order (NLO) pQCD calculations. The present analysis follows 
a data-driven approach with minimal theoretical assumptions where, in order to minimize systematic un- 
certainties, only measurements by ALICE are considered. In this approach, the pp reference spectrum is 
obtained by interpolating the differential yields d^Nj!^ /drjdpT of charged particles measured in inelastic 
pp collisions at y/s = 0.9 and 7 TeV by ALICE |[T9ll20l . The interpolation is performed in bins of pj, 
based on the assumption that the increase of the yield with y/s follows a power law. Above pr =2 GeV/c, 
the measured spectra at the two energies are parametrized by a modified Hagedorn function 1(211 and a 
power law to reduce bin-by-bin fluctuations. Systematic uncertainties on the pp reference spectrum arise 
from the experimental errors of the measured spectra at 0.9 and 7 TeV, from the parametrization, and 
from the interpolation procedure in y/s. The combined statistical and systematic data enws result in a 

9- 10% uncertainty on the pp reference spectrum at y/s = 2.76 TeV, depending on pr- The interpolation 
procedure was verified using PHOJET 121 and PYTHIA [23] (tunes D6T El and PerugiaO |[25l) at 0.9, 
2.76 and 7 TeV. The generated and interpolated spectra at 2.76 TeV agree within the quoted uncertainties. 
Finally, the scaled pp yield in a given centrality class is obtained by multiplication of the pp reference 
spectrum with (A'coii), see Table [T] The uncertainty in (A'^coii) results in an additional /^j- -independent 
scaling uncertainty on the scaled pp reference. 

Alternative approaches to derive the pp reference spectrum are investigated to study the sensitivity of Raa 
to the specific choice of our method. Replacing in the interpolation the pr spectrum at 0.9 TeV by the 
one measured in pp at y/s = 1 .96 TeV in | T] | < 1 by the CDF collaboration |26il results in a pp reference 
spectrum which is 5-15% lower than the reference spectrum described above. A different procedure to 
obtain a pp reference is based on a scaling of the pj spectra at 0.9 or 7 TeV to 2.76 TeV by the relative 
dependence predicted by NLO pQCD calculations ETl (referred to as "NLO scaling"). Using the 
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Figure 2: The pj distributions of primary charged particles at mid-rapidity {\ri\ < 0.8) in central (0-5%) and 
peripheral (70-80%) Pb-Pb collisions at ^/s^ = 2.76 TeV. Error bars are statistical only. The systematic data 
errors are smaller than the symbols. The scaled pp references are shown as the two curves, the upper for 0-5% 
centrality and the lower for 70-80%. The systematic uncertainties of the pp reference spectra are contained within 
the thickness of the line. 

7 TeV spectrum as a stalling point, good agreement with the reference obtained from interpolation is 
found. Starting instead from 0.9 TeV results in a spectrum which is 30-50% higher than the interpolation 
reference. The pp reference spectra derived from the use of the CDF data in the interpolation and from 
NLO scaling of the 0.9 TeV data are used in the following to illustrate the dependence of Raa at high pr 
on the choice of the reference spectrum. 

The pt distributions of primary charged particles in central and peripheral Pb-Pb collisions at 2.76 TeV 
are shown in Fig. [2l together with the binary-scaled yields from pp collisions. The -dependence is 
similar for the pp reference and for peripheral Pb-Pb collisions, exhibiting a power law behaviour at 
Pt > 3 GeV/c, which is characteristic of perturbative parton scattering and vacuum fragmentation. In 
contrast, the spectral shape in central collisions clearly deviates from the scaled pp reference and is closer 
to an exponential in the pj range below 5 GeV/c. 

Figure |3] shows the nuclear modification factor Raa for central and peripheral Pb-Pb collisions. The 
nuclear- modification factor deviates from one in both samples. At high pr, where production from hard 
processes is expected to dominate, there is a marked difference between peripheral and central events. In 
peripheral collisions, the nuclear modification factor reaches about 0.7 and shows no pronounced pj de- 
pendence for Pt >2 GeV/c. In central collisions, Raa is again significantly different from one, reaching 
a minimum of Raa ~ 0.14 at pT = 6-7 GeV/c. In the intermediate region there is a strong dependence 
on Pt with a maximum at pT = 2 GeV/c. This may reflect a variation of the particle composition in 
heavy-ion collisions with respect to pp, as observed at RHIC |[28ll29l . A significant rise of Raa by about 
a factor of two is observed for 7 < < 20 GeV/c. Shown as histograms in Fig.[3l for central events only, 
are the results for Raa at high pt, using alternative procedures for the computation of the pp reference, 
as described above. For such scenarios, the overall value for Raa is shifted, but a significant increase of 
Raa in central collisions for pt >^ GeV/c persists. 

In Fig. |4] the ALICE result in central Pb-Pb collisions at the LHC is compared to measurements of 
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Figure 3: Raa in central (0-5%) and peripheral (70-80%) Pb-Pb collisions at ^/s^ = 2.76 TeV. Error bars 
indicate the statistical uncertainties. The boxes contain the systematic errors in the data and the pr dependent 
systematic errors on the pp reference, added in quadrature. The histograms indicate, for central collisions only, 
the result for Raa at pr > 6.5 GeV/c using alternative pp references obtained by the use of the pp measurement 
" = 1 .96 TeV ||26| in the interpolation procedure (solid) and by applying NLO scaling to the pp data at 0.9 
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Figure 4: Comparison of Raa in central Pb-Pb collisions at LHC to measurements at y/s^ ~ 200 GeV by the 
PHENIX ilOl and STAR ^ experiments at RHIC. The error representation of the ALICE data is as in Fig. [3] 
The statistical and systematic errors of the PHENIX data are shown as error bars and boxes, respectively. The 
statistical and systematic errors of the STAR data are combined and shown as boxes. The vertical bars around 
Raa = 1 indicate the pr independent scaling errors on Raa- 
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Raa of charged hadrons = 200 GeV) by the PHENIX and STAR experiments |[30l[3Tl at RHIC. 

At 1 GeV/c the measured value of Raa is similar to those from RHIC. The position and shape of the 
maximum at 757 ~ 2 GeV/c and the subsequent decrease are similar at RHIC and LHC, contrary to 
expectations from a recombination model fyT\. Despite the much flatter pj spectrum in pp at the LHC, 
the nuclear- modification factor at pj = 6-7 GeV/c is smaller than at RHIC. This suggests an enhanced 
energy loss at LHC and therefore a denser medium. A quantitative determination of the energy loss 
and medium density will require further investigation of gluon shadowing and saturation in the present 
energy range and detailed theoretical modeling. 

In summary, we have measured the primary chai^ged particle pj spectra and nuclear modification factors 
Raa in central (0-5%) and peripheral (70-80%) Pb-Pb collisions at y/i;;;^ = 2.76 TeV with the ALICE 
experiment. The nuclear modification factor in peripheral collisions is large and independent of pj for 
Pt GeV/c, indicating only weak parton energy loss. For central collisions, the value for Raa is found 
to be ~0.14 at pj = 6-7 GeV/c, which is smaller than at lower energies, despite the much less steeply 
falUng pt spectrum at the LHC. Above 7 GeV/c, Raa increases significantly. The observed suppression 
of high Pt particles provides evidence for strong parton energy loss and large medium density at the 
LHC. 
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